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Abstract: We generalize and simplify an earlier approach. In three dimensions we 
present the most general averaging formula in lowest order which respects the require- 
ments of covariance. It involves a bitensor, made up of a basis of six tensors, and contains 
three arbitrary functions, which are only restricted by their behavior near the origin. The 
averaging formula can also be applied to the Einstein tensor. If one of the functions is 
put equal to zero one has the pleasant property that the Einstein tensor of the averaged 
metric is identical to the averaged Einstein tensor. We also present a simple covariant 
extension to static perturbations in four dimensions. Unfortunately the result for the 
Einstein tensor cannot be extended to the four dimensional case. 
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1 Introduction 



The energy-momentum-tensor used in cosmological models is an average over the non 
homogeneous tensor present in nature. Therefore one needs an averaging prescription for 
the energy momentum tensor and for the metric. This provides a fundamental problem 
because of the freedom of choice of coordinates. The averaging problem in general 
relativity was first raised by Shirkov and Fisher [1] in 1963. The authors of [1] suggested 
to integrate the metric tensor over a four dimensional volume with the familiar factor 



\J—g in the measure. Such an expression is, however, not covariant due to the freedom of 

performing local transformations. A covariant averaging prescription can be constructed 
by introducing a bivector g@(x, x') of geodesic parallel displacement, as discussed in the 
appendix of [2] . This transforms as a vector with respect to coordinate transformations 
at either x or x' and maps a vector Ap{x') to A a {x) = gP(x,x')Ap(x'), analogously 
for higher order tensors. An averaging with the help of bivectors was also used in the 
work of Zalaletdinov [3] where the emphasis was on the commutativity of averaging and 
covariant differentiation. As remarked by Stoeger, Helmi, and Torres [3], the method 
of using a covariantly conserved bivector is not applicable to the metric, because the 
covariant derivative of the metric vanishes. The metric is therefore invariant under this 
averaging procedure. In the thesis of Behrend [5] the metric is represented by tetrads 
and the averaging performed over the latter. The tedrads are chosen according to a 
covariant minimalization prescription. 

This is only a very brief survey of the literature. For more references, as well as 
the implications for the fitting problem, back reaction, contributions to dark energy, we 
refer e.g. to the monograph of Krasinski [6] and to the comprehensive recent reviews of 
Buchert [7] and of Malik and Wands 0. 

In a previous paper [H] we gave an explicit solution. The derivation was rather com- 
plicated, furthermore we could not classify the complete set of solutions. The present 
approach is much more transparent and straight forward. It naturally leads to a classi- 
fication of all solutions which fulfill the central requirements (1.1) and (1.2) below. 

Under a covariant averaging process we understand a prescription which has the 
following properties. Let two observers describe the same physics in different coordinate 
systems S and 5", with metric tensors g^ v and g' . Both of them apply a definite 
averaging procedure in their respective systems, resulting in the averaged metrics 
< gp,v > an d < g' >, respectively. Then the results have to be connected by the same 
transformation as the original metric, i. e. 



In other words, the operations of averaging and of coordinate transformations have to 
commute. 

Furthermore, averaging over a region which is closely located around some point 
should, of course, reproduce the metric at this point. In this situation only the metric 
at the origin is relevant. This means that a constant metric has to be reproduced by the 
averaging process, i. e. 




< g'av >=< 5V >' ■ 




< Vtiv >= V/iu for constant rj^. 
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It is rather obvious that space and time cannot be treated at the same footing in an 
averaging prescription. Therefore, as usual, we have to assume that a reasonable foliation 
into space and time can be performed. 

In sect. 2 we present the most general lowest order three dimensional covariant 
averaging formula. The mapping of the metric gki{x!) to the averaged metric < g mn > (x) 
is represented by a bitensor K^Jjs! — x) which we will specify in detail. A product of 
bivectors, as frequently used in the literature, is not sufficient. The bitensor K® n (x! — x) 
is a superposition of a basis of six tensors, and contains three arbitrary functions u(r), 
v(r), and w(r), which depend upon the distance r = |x' — x| and are only restricted by 
a prescribed behavior near the origin. They are necessarily singular there. Our earlier 
formula in [9] is a special case of our general formula as it should be. In sect. 3 we 
show that, under the condition w(r) = 0, we can use the same formula which works for 
the metric as well for the Einstein tensor. This implies that the Einstein equations for 
the averaged metric are identical to the averaged equations. In sect. 4 we show why 
iteration of the averaging formula does not make sense. Therefore one has to choose 
reasonable functions from the beginning. In sect. 5 we discuss a simple four dimensional 
covariant generalization for static perturbations. The proof for the covariant averaging 
of the Einstein tensor "almost" goes through also in this case, but fails at the very end. 

We claim that our formulae present the most general framework for a lowest order 
averaging prescription which respects (1.1) and (1.2). 

2 The most general lowest order three dimensional 
covariant averaging formula 

Consider a perturbed flat metric «7fcz(x) = 5ki + fow(x). In lowest order there must be 
a linear relation between the original metric and the averaged one. For simplicity, we 
consider the averaging at the origin for the moment. 

<<7»»>(0) = /ii^(x)s w (x)^. (2.1) 

The same relation holds for the perturbation h mn , because, as we shall see, averaging 
of the constant S mn gives back this constant as required by (1.2) (the above K™ n (x) is 
identical to the lf*J n (x) in [9]). 

In order to guarantee covariance with respect to rigid rotations, one has to construct 
the most general structure of the tensor (tensor in the sense of linear algebra) K^ n {x) 
which is symmetric under the exchange m «-> n and under k <-> I. This is represented by a 
complete set of six independent tensors which are multiplied by functions A(r), • • • , F(r), 
which depend only upon the distance r = |x|: 
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.(2.2) 



We emphasize that this is the most general tensor structure which one can write down. 
Naive averaging would only make use of A{r), normalized to J °° 2A(r)r 2 dr = 1, while 
all the other functions B(r), ■ ■ ■ , F(r) vanish. 

Let us now apply an arbitrary infinitesimal transformation x k = x' k + C, k , which leads 
to a change dg^i = £k,i + £i,k ~^ %£k,h when contracted with the symmetrical tensors in 
(2.2). Covariance according to the requirement (1.1) implies that the change of the rhs 
of (2.1) must be identical to the change of the lhs, i.e. to £ m , n (0) +£ n ,m(0). This means 
that the rhs cannot depend upon £fc(x), except at the origin x = 0. To see how this can 
happen we perform a partial integration with respect to xj. The change of the integrand 
on the rhs of (2.1) then becomes —2K™ n ,i £ fc . 

This has to vanish for arbitrary which implies K™ n ,i = 0. Therefore K™ n has to 
be a curl with respect to I, and, by symmetry, also with respect to k, i.e. we can put 



The tensor T^ d n (x) has a decomposition analogous to (2.2), we denote the six radial 
functions by the corresponding small letters a(r), • • • , f(r). 



Tl(x) = a(r)M + M+^)^ M + c(r)5 mn ^ + rf(r)^^ 

ry ry d ry ry^ ry ry ^ ry ry ^ ry ry ry^ ryd 

+e(r)[5 m — + 5 m — + K^r + 5 n—\ + f( r ) ~ 4 ■ ( 2 - 4 ) 



This can be introduced into (2.3) and, after performing the differentiations, be compared 
with the general decomposition (2.2). The result is 




(x) = 44^ c T^(x). 



(2.3) 



D(r 



B(r 



E(r 



C{r 



A(r 
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Up to now we have six functions a(r), • • • , f(r). But only three of them are relevant. To 
see this we first put a'(r) = a(r)/r and b'(r) = b(r)/r. One then finds that e(r) and a(r) 
only appear in the combination e(r) — 5(r)/2, while c(r) and 6(r) only appear in the 
combination c(r) — b(r). One thus can rename these combinations as e(r) and c(r), i. e. 
one can put a(r) = b(r) = in (2.5). The next step is more subtle. Introduce two new 
functions g(r) and w(r) by g(r) = [c(r) +/(r)]/2 and w(r)/r 2 = C(r) — D(r). Eliminate 
c(r) = [-rV(r) /4 - rd'(r) /4 + d(r)] + [rg'(r) /2 + #(r)] - iu(r) /4, and 
/(r) = [r 2 d"(r) /4 + rtf (r) /4 - d(r)] - [rg'(r) /2 - g(r)) + w(r) /4, 
and use the functions g(r) and w(r) instead of c(r) and /(r). If we further put 
g(r) = rg'{r) — 2g(r), and eliminate d(r) by introducing the function u(r) = d(r) — 2g(r), 
we find that e(r) and g(r) only appear in the combination v(r) = e(r) —g(r). Therefore 
we end up with three relevant functions u(r),v(r),w(r). (These u(r),v(r) have, of course, 
nothing to do with the projective coordinates u, v introduced in [9].) We thus have found 
the following representations for the functions A(r), ■ ■ ■ , F(r) in (2.2): 



Afy) = 4KM + ^] + 2[^-^]4^, 

B(r) = i [rB >) +5u >)-^]-4[^-^l + i[^) + 2^1, 

C(r) = - i K» + 3 W "(r)-3^]+4[^-2^]- i ^, 

D( r ) = __[ rM '"( r ) + 3 M "( r )_3^_^] +4 [^_i_2^]--[^^ + 4^| 



C(r 



u>(r) 



r 
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E(r) = i[ M >)-^]-2[^-2^] + i4 :1 



2 r r r 2 2 r 2 
. .., --^-16-^ +T—^- 4 -i 



F(r) = 7 [™>)-5 M ''(r) + 13^-16^] + ^[^-4^]. (2.6) 



We next demand that the averaging over a region which is closely localized around 
the origin should give back the metric g mn {0) at the origin. In this case we may put 
<7fcj(x) = <7m(0) in (2.1) and take it out in front of the integral. The angular averages 
J dQ/in can be performed using 

^-~<W ^^^^ m J kl + 6 k m 6 l n + 6 k JJ etc. (2.7) 
r z 3 r 4 15 

Comparing the factors in front of the terms <? mn (0) and g\{Qi)5 mn on both sides one thus 
obtains 



r°° 4 2 f 00 1 1 1 

1 = / \2A{r) + -£(r) + —F(r)}r 2 dr, = / [B(r) + -C(r) + -D(r) + —F(r)]r 2 dr. 
Jo o 15 Jo 6 6 15 

(2.8) 

This derivation also shows the validity of (1.2), i. e. that the averaging of a constant 
gives back this constant. Therefore our averaging formula may be applied to the metric 
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g mn , as well as to the perturbation h mn . The relations can also be rephrased in the 
property 

fK^)^ = \[S k J l n + S k n S l J. (2.9) 

Inserting the representations (2.6) one finds that the integrands in (2.8) can be written 
as derivatives. Assuming that there are no boundary terms at infinity one thus obtains 



1 = [-^rV(r) + -L r V(r) + ^rw(r) - ^rv(r) - ^ra(r)] r=0 , 

1 11 4 4 1 

= [- — rV(r) - — rV(r) + —ru(r) + -rv(r) - — rw(r)] r=0 . (2.10) 

1U oU 10 o 1U 

The derivation above makes clear that the representation (2.1), (2.2), together with the 
form (2.6) and the boundary conditions (2.10), is the most general first order covariant 
averaging formula in three dimensions. 

Before we discuss the implications of (2.10) we consider the boundary terms at the 
origin which arise from the partial integration of 2£k,i, where was the shift of an 
infinitesimal transformation. Restricting the integration to the outside of a small sphere 
of radius e these are just the surface terms which arise from GauB's theorem. Using 

*^(x)&,i (x)^ = -2 / iOxA^. (2.11) 
Air J r =e e Air 

Because invariance with respect to translations and to rigid rotations around the point 

of consideration is manifest, one can restrict to transformations which leave the origin 

fixed, such that £fc(x) = ^(O)^ + 0(e 2 ). If we insert this into (2.11) and use the 

representation (2.2) for K™ n (x) we can perform the angular averaging. The result has 

to be identical to the change £ m , n (0) + £ n ,m(0) of the lhs. Comparing the factors of 

£m,n(0) +£n,m(0) and of 25 mn ^(0) on both sides one obtains 



1 = 
= 



2r 3 

15 



[5A(r) + D(r) + 7E(r) + F(r)} 



r=0 



— [5B{r) + 5C(r) + D(r) + 2E(r) + F(r)} 
15 



r=0 



(2.12) 



If one introduces the representations (2.6) for A(r), • • • , F(r), one obtains again the 
conditions (2.10). 

The boundary conditions (2.10) have drastic consequences for the behavior of the 
functions at the origin. They imply 



u(r) = U[-i]/r + 0(l),v(r) = + 0(1), w(r) = W[-i]/r + 0(1), with 

5 1 13 1 



u l-i] = 4 + g w l-i] 



(2.13) 
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This result is unpleasant but unavoidable. The functions u(r), v(r), w(r) are singular 
and behave like 1/r near the origin. This implies that the functions A(r), • • • , F(r) go 
like 1/r 3 . At first sight this might look as if the integrand in (2.1) is not integrable at 
the origin. This is, however, not the case. Expand ^/(x) = gki(0) +0(x). The constant 
term gki(0) can be taken out of the integral, the angular averaging of the l/r 3 -term 
vanishes. This constant term has just been treated in detail. The rest is of order r/r 3 
and therefore integrable in three dimensions. 

Nevertheless one would have preferred functions A(r), • ■ ■ , F(r) in the averaging for- 
mula which are smooth at the origin. But covariance definitely prohibits such a smooth 
behavior. If one considers the change of the integral which arises from an infinitesimal 
coordinate transformation x k = x' k + the boundary terms of the partial integration 
have to reproduce the change of the metric at the origin, i.e. £ m , n (0) + £ n ,m(0). This 
enforces the singular behavior of the functions. Smooth functions could not produce 
boundary terms. 

In [H] we presented a special solution for a covariant averaging procedure. It con- 
tained a function f(r) (which has nothing to do with the /(r) in (2.4)), normalized 
to Jq° f(r)dr = 1, as well as two integrals of f(r), namely F(r) = — / r °° f{r')dr' and 
G(r) = — J r °° f(r')/r'dr'. This solution must be a special case of our general formula. To 
demonstrate this we have to introduce another integral H(r) = — J r °° f(r')/r' 3 dr'. Then 
we obtain our old solution if we put 



«(r) = + f G(r) - B ~r>H(r), v(r) = H£M _ ^H(r). (2.14) 

When making this comparison one has to take take care of the correct factors r 2 . In 
[9] we used the integration element drdQ/Air because we had to perform several partial 
integrations with respect to r there, while in (2.1) we use d 3 x/4n = r 2 drdfl/4:7r. 



3 Covariant averaging of the Einstein tensor 

Besides the metric, the Einstein tensor is the most important object in general relativity 
because it enters, together with the energy momentum tensor, directly the field equa- 
tions. It would be highly desirable if one could average the Einstein tensor in exactly the 
same way as the metric tensor, and if the averaged Einstein tensor would be identical to 
the Einstein tensor derived from the averaged metric. The problem that the averaged 
equations are not identical with the equations for the averaged metric would then dis- 
appear. We repeat and extend some of the steps of [9] in order to make the paper self 
contained. 

In first order of the perturbation the Einstein tensor becomes 

on — h 1 -\-h i —h i —h i —h ij r) -\- h ■ ■ ij 1 
ilJ mn 1 L iimn\ 1 L mnii IL irm n IL ini m u mn \ 'Hjt u mn- \°* / 

Indices are raised and lowered with 5{j here, so their position is in fact irrelevant. The 
averaging formula (2.1) is now used for an arbitrary point x, the integration variables 
are denoted by a prime, and the tensor K^ n (x. f — x) depends on the difference x' — x. 
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The distance r' now means r' = |x' — x|. Due to the singular behavior of K!^ n (x.' — x) 
at zero distance one has to treat this region separately. Let us assume that we calculate 
the Einstein tensor near the origin, and expand 

M*') = Mo) + K, a (OK + \Ki? h WA + M*')- (3.2) 

The term ^jj(x') is of order x' 3 , no problems with potentially divergent contributions or 
boundary terms from partial integrations can arise. We begin with this term. 

The first possibility is to average the metric in the way described before. Subsequently 
one calculates the Einstein tensor from (3.1), using the averaged metric on the rhs. The 
factor of h k i(x.') in the integrand then becomes 

The second possibility is to calculate the Einstein tensor Gki with the old metric and 
then average it with our formulae in exactly the same way as we averaged the metric 
tensor. Shift the partial derivatives from the perturbation h to K, and rename dummy 
indices where necessary such that h k i appears in all six terms. The factor of /ifcz(x') m 
the integrand now becomes 

jkl Tsij ski _|_ jskl i l(ki I jy'li k jy-i j rfe/ , jy-i kl /n A \ 

a ran mniij ' mnii mnii mnii mniij ^ ' mnii " / 

Consider the difference A mnk i = I mn ki — Jmnki- If one combines the terms appropriately 
one has 



A-mnkl (-K ' ikliirin ^mniikl ) {-^imklin ^mnkiil ) {J^inkli m ^mnliik ) 

-( K lkh] Smn - Sia) + {K ijkh lJ 5 mn - K mnij , lJ 5 H ). (3.5) 

The further investigation can be greatly simplified if one decomposes 

K mn kl = K\r} n kl + ^Llkli (3-6) 

where K^} nkl is symmetric against the exchange (m, n) <-> (k, I) and K^ nkl antisymmet- 
ric. Obviously K^} nkl consists of the terms with A(r), B(r), E(r), F(r) in (2.2), together 



[A] 

mnkl 



with the symmetric combination [C(r) + D(r)][5 mn x k xi/r 2 + 5kix m x n /r 2 }/2, while K 
consists of the antisymmetric combination [C(r) — D(r)][5 mn XkXi/r 2 — 5kix m x n /r 2 }/2. 
From (3.5) it is seen that the symmetry relations in A mnk i are just reversed with respect 
to K mn ki, i.e. the symmetric part A^ nkl is obtained from K^l kl , while the antisymmet- 
ric part A^]^ is obtained from K^} nkl . Clearly A mnk i has a decomposition analogous to 
K mn ku with coefficients A(r), ■ ■ ■ , F(r), say. 

Let us first investigate the antisymmetric part A^ nkl which arises from K^} nkl . Know- 
ing the structure A^ nkl = [C(r) - D(r)\ [5 mn x k xi/r 2 - 5 H x m x n /r 2 }/2, we can simplify 
the investigation by taking the trace k — I, thus we only need to calculate 
= [C(r) - D(r)][5 mn - 3x m x n /r 2 ]/2. The result is 
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C(r) - D(r) 
2 ~ 

-[B"{r) + C "( r ) + D "( r ) ] + [s / (r ) + F / (r) ]I + [c(r) + £, (r) + 2F(r)]i 
2 r v 

= w"{r)/2r 2 -u/(r)/r 3 -2w(r)/r 4 . (3.7) 

The functions w(r) and v (r) have dropped out completely, the expression vanishes if 
w(r) = 0. To see that this condition is also necessary we only need to consider, e.g. the 
contractions m = n, k = I of the symmetric part A m ^ k i, which gives 

A£[ mfc = -2w"(r)/r 2 - 2w'(r)/r 3 + 8w(r)/r\ (3.8) 

Both (3.7) and (3.8) have to vanish and this is the case if and only if w(r) = 0. This is 
equivalent to the equation C(r) = D(r) } i.e. to the symmetry relation K^ nkl = 0. 

We have to check that the independence of the order of averaging also holds for 
the first three terms in the decomposition (3.2). For the constant and linear terms 
hij(0)+hij, a (0)x' a this is trivial, they do not contribute to the Einstein tensor, irrespective 
of the order of averaging. In the average of the quadratic term hij, ab (0)x' a x' b /2, we 
substitute x' — x = y, such that x' a x' b = x a X\, + x a yb + Da%b + VaVb- In the first term 
one can take out x a x b in front of the integral, the remaining integral is given by (2.9), 
therefore hij, ab (0)x a Xf ) /2 is reproduced. The other three terms which are linear and 
constant with respect to x do not contribute to G mn . If, alternatively, we first calculate 
Gran-, which is constant for this contribution, it is as well preserved by the averaging. So 
we have seen by direct evaluation that the possibly dangerous low order terms in (3.2) 
do not generate problems. 

We finally have shown that the Einstein tensor of the averaged metric is identical 
to the averaged Einstein tensor of the original metric if and only if wir) = 0, which is 
equivalent to C{r) = D(r), i.e. the symmetry relation K mn ki = Kki mn - In [9] we had 
assumed this symmetry in order to simplify the discussion, we now have shown that this 
condition is necessary. In the following we will always assume w(r) = 0. 

An important feature for the understanding of this property are the symmetry rela- 
tions shared by our averaging formula and by the Einstein tensor. This becomes clear if 
one writes 

2G mn = T mn h k i 1 (3-9) 

with the operator 



Tl = 5 kl d m d n + -(6^6 l n + 6 k Jl)m--(6td n d l + 6ld^ 

-5 kl 5 mn d l d % + 5 mn d k d l . (3.10) 

Both tensors, K mnk i as well as T mn M, are symmetric under m <-> n, under k <-> I, and 
under (m,n) (k,l). These symmetries were essential in order to show the vanishing 
of the difference (3.5). 
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The result for the covariant averaging of the Einstein tensor is certainly not trivial. 
For the Ricci tensor, which does not fulfill the above symmetry properties, the relation 
is not valid. 

4 Iteration and stability 

The general behavior of iterations is most easily first studied in a simple one dimensional 
toy model. Consider an averaging formula 

<g>(x) = 9[1] (x) = j ' K{x - x')g{x')dx' ', (4.1) 

with K(x) real and even, and normalized to / K(x)dx = 1. Because (4.1) is a convo- 
lution, it is convenient to work with the Fourier transforms g(p) — j g(x)e~ tpx dx, etc. 
This implies K(0) = 1 and g\i]{p) = K(p)g(p). The iteration of order n of the averaging 
procedure becomes 

g [n] (p) = K n (p)g(p). (4.2) 

From this it is immediately clear that iteration does not make much sense. If > 1 

for some values of p, the iteration will diverge. If < 1 the iteration will converge 

to for these p. A stable averaging prescription K{p) 2 = K{p) will be obtained if and 
only if K(p) only takes the values or I. Thus consider a (finite or infinite) sequence 
do = < b < a\ < b\ < • • •, and put K{p) = Xi„6(a n < \p\ < b n ). This implies the 
general form K(x) = (1/irx) X] n ( sm ^n^ — sina„a;) for a stable averaging function in this 
simple toy model. 

One can start a similar investigation for our three dimensional averaging formula. 
It is easy to formulate the iterations in Fourier space, but it appears hard, probably 
impossible, to fulfill the conditions for a stable solution together with the representation 
(2.6) and the boundary conditions (2.13). The requirement of covariance prevents a 
stable averaging procedure. 

5 Static perturbations in Minkowski space 

Our extension to the four dimensional case is rather modest but practical. We consider 
a Robertson Walker metric with k — 0. A substitution r = r'/a(t), with a(t) the cosmic 
scale factor, brings the line element into the form ds 2 = dr 2 — dt 2 + • • •, where the 
corrections are small as long as the region of averaging is small compared to the Hubble 
length. Therefore we can use the Minkowski metric (1,1,1,-1) as the unperturbed metric. 
We further assume that the perturbation is approximately static. To keep this situation, 
we restrict the admissible transformations to rigid translations, rigid spatial rotations, 
and infinitesimal transformations which keep the time unchanged. This means that 
£° = 0, and £ m is independent of t. Furthermore we can drop all time derivatives in the 
perturbed metric. Under these restrictions the perturbations hoo and h m o become gauge 
invariant. 

We average the perturbation with the following simple ansatz. 
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</w>(x) = jifKx'-^^x'jp (5.1) 

</^oo>(x) = /P(x'-x)/ i00 (x')^, (5.2) 

Q(x'-x)/i m0 (x') — . (5.3) 

Here P(x' — x) and Q(x.' — x) have to be rotation invariant and correctly normalized. 
Equations (5.1) - (5.3) are the simplest generalization of our previous formula. For static 
perturbations they are covariant in the sense of (1.1) with respect to static transforma- 
tions. 

Let us now consider the (lowest order) Einstein tensor which, in the static case, reads 



9(7 — oni s ) 4- h° —h° ' l S C5 A) 

^^mn — ran ' ' h Q-imn "'OH u mnt \°-^) 

2Goo = hii-hij? (5.5) 
2G m0 = h m0 ,l —h i0 , l m . (5.6) 

Here G$ n is the spatial part of the Einstein tensor in (3.1). Under the assumptions above, 
the additional terms in G mn , as well as G o and G m0 are invariant under infinitesimal 
static transformations. 

We investigate whether the Einstein tensor of the averaged metric can be identical to 
the averaged Einstein tensor. Although everything looks promising at the beginning, we 
will obtain a negative answer at the end. One may therefore skip the rest of this section. 

As in the previous section we treat the constant, linear, and quadratic terms which 
were split off in (3.2) separately. Again for these the averaging is independent of the 
order in which it is performed. We start with G mn . For G$ n we have shown in the 
previous section that the result is independent of the order of averaging. We can restrict 
to the additional contributions in (5.4). 

If we average the perturbation h® according to (5.2) and introduce into (5.4) we obtain 
the integrand [P, mn — P,\ S mn ]hQ. If, alternatively, we first calculate 2G mn in (5.4) with 
the old metric and then average it in the same way as (5.1), i.e. replace hki by G^ there 
on the rhs, and shift the partial derivatives from the perturbation to the multiplying 
functions, we obtain the integrand [K^ n:k i — 

This gives the condition 

p — P*A — K kl ,,—K kl (^7\ 
1 imn 1 ii u mn lv mmkl sy mnktl • \' J -') 

We next apply the same procedure to Goo- If we average the perturbations h\ and hij 
according to (5.1) and introduce into (5.5) we obtain the integrand \K %kl ^ —K k J, l ^]hki- 
If, alternatively, we first calculate 2G o i n (5-5) with the old metric and then average it in 
the same way as (5.2), i.e. replace /i o by Goo there, and shift the partial derivatives from 
the perturbation to the multiplying functions, we obtain the integrand [P,j S kl — P, kl }hki- 
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The condition K\ kl ^ —K^j,^ = P, J j 5 kl — P, kl which now arises is identical to (5.7) if 
one renames the dummy indices k, I, subsequently replaces m,n by k, I, and uses the 
symmetry K mn ki = Kki mn - We thus only need to consider (5.7) in the following. 

An elementary calculation under proper consideration of the singular behavior at the 
origin gives (of course we demand C(r) = D(r), i.e. w(r) = in this section) 



KLm -K k mnk !i = ^i(r)= + K 2 (r)5 mn - -[d rn d n - 5 mn A]5^(r), (5.8) 



with 



+2 |^i) + 4 £W, 

K2(r) = V(r) + ^l-2[B-(r) + 2^] + 2[^i-^i] + 4^M 

ry /y /y ry A ry 

(5 . 9) 

The term with the 5-function is most conveniently obtained by multiplying (5.8) with x l x^ 
and integrating. If one introduces the representation (2.6) and considers the behavior 
(2.13) near r = 0, one finds that the singular terms in u(r) and v(r) cancel, Ki(r) 
and K2(r) behave at most like 1/r 4 for small r. By a suitable choice of the remaining 
freedom in u(r) and v(r) one could also remove any singularities, but the convergence of 
the integrals used in the following is guaranteed anyhow. If one equates the factors of 
-[x m x n /r 2 + 5 mn ]/r, of [x m x n /r 2 - 5 mn ] , and of [d m d n -<5 mn A]5 (3) (r) in (5.7) one obtains 



with 



P(x) = P(r)-y5 (3) (r), (5.10) 



P'(r) = - r -[ Kl (r) + K 2 {r)l P"(r) = ±[ Kl (r) - K 2 (r)]. (5.11) 
This implies the integrability condition 

K[(r) + K>(r) + ^)=0. (5.12) 

r 

It is fulfilled if one introduces the representations (5.9) and (2.6). The normalization of 
P becomes 



roc „ 1 poo „ 1 poo 3 

/ r 2 P(r)dr = -- / r z P'(r)dr = - / r^K^r) + K 2 (r)]dr = -. (5.13) 
Jo 3 jo 6 Jo 2 

From (5.10) this implies the correct normalization of P. The last integrand in (5.13) (as 
well as the integrand with K\(r) or K 2 (r) alone) turns out to be a derivative when one 
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introduces (5.9) and (2.6). Therefore no freedom is left, the integral is determined by 
the boundary term at zero, and the latter is fixed by the behavior of u(r) and v(r) in 
(2-13). 

This result is unpleasant. It fixes the integrals over r A Ki(r) and r i K 2 {r) ) and leads 
to the unwanted 5-function contribution in (5.8) and (5.10). The 5-function in (5.10) 
would imply that < h$> (x) would contain a contribution — }iq(x)/2, i. e. a contribution 
which is not averaged, and furthermore, with opposite sign. Such a contribution cannot 
be tolerated. 

One could use a more general ansatz in (5.1), (5.2), where terms with Skiti® and 
(xkXi/r 2 )tiQ are inserted into the rhs of (5.1), and a gauge invariant combination of ft| 
and (x k x l /r 2 )h k[ into (5.2). We found that this does not help to solve the problem. At 
the end the relevant integrands again turn out to be derivatives, and everything is fixed 
by the boundary conditions. We conclude that the previous result about the covariant 
averaging of the Einstein tensor cannot be extended to four dimensions. 

Therefore we keep things simple and work with (5.1) - (5.3). We may choose rather 
arbitrary functions P and Q in (5.2), (5.3), of course with the correct normalization and 
somehow related to the functions in (2.2). Nevertheless we favor the choice w(r) = 
also in the four dimensional case. 

6 Outlook and conclusions 

In this paper we presented a covariant averaging prescription which fulfills two essential 
requirements. 

• Under coordinate transformations the averaged metrics are connected by the same 
transformation as the original metrics, i.e. < g' >=< g^ v >' (eq. (1.1)). 

• The averaging of a constant metric reproduces this metric, i.e. < r]^ >= rj^ for 
constant rj^ (eq. (1.2)). 

In three dimensions we gave the complete solution of the problem. There is a linear 
connection between the original perturbation /ijw(x') and the averaged < h mn > (x), 
represented by a tensor (tensor in the sense of linear algebra) K® n (x! — x). (Such 
a connection has also been discussed by Boersma [10J, although without going into 
details.) The tensor K^ n (pc' — x) is a superposition of a basis of six bitensors which are 
symmetric with respect tom«n and to k I. A product of bivectors, as sometimes 
suggested in the literature, is not sufficient. The representation contains three functions 
u(r),v(r),w(r) which depend upon the distance r = |x' — x|. They have to be singular 
in a definite way at the origin, in order to fulfill the requirement (1.1) of covariance. For 
w(r) = we found the welcome property that the Einstein tensor can be averaged in the 
same way as the metric, and that the Einstein tensor belonging to the averaged metric 
is identical to the averaged Einstein tensor. 

We have further seen that it does not make sense to iterate the averaging procedure. 
Therefore one has to choose a reasonable ansatz for the functions u(r) and v(r) from the 
beginning, of course with the correct boundary conditions. Three simple choices suggest 
themselves: 
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Exponential function: 



5 e~ r/r ° 13 e~ r/r ° 

"'esepv ) 

Gaussian: 



u exp (r) = l-^—i v exp (r) = -gg— — ■ (6-1) 



5 e -(r/r ) 2 ^ ^ 13e~^ //r °' ) 

Averaging over a sphere: 



5(r -r) 4 13 (r -r) 2 

U S phere{r) = ~ J 9(r - r), tVfterel 7 *) = ~™ 2 ( r O ~~ ( 6 " 3 ) 

In the last case we took care not to get 5— like contributions in the derivatives from the 
boundary at r = r . 

The functions A(r), • • • , F(r) can be easily obtained from this with the help of (2.6), 
there is no need to show the explicit expressions here. 

The generalization to static perturbations in Minkowski space, for many applications 
an excellent approximation to the realistic case, is quite simple. The property for the 
covariant averaging of the Einstein tensor could, however, not be generalized to the four 
dimensional case. 

One can hope to obtain a sufficiently smooth metric after performing the average. 
But it is important to note that we are still free to perform gauge transformations. By 
an unfavorable choice of gauge the "smoothed" metric can become wavy and irregular. 
All one can achieve is that the final metric becomes equivalent to a smooth metric. 

The present investigation was already quite elaborate, it gave the mathematical 
framework for a covariant averaging prescription. Applications have to be postponed 
to forthcoming work. The central question to be investigated is, of course, how far the 
averaging of inhomogeneities can mimic the presence of dark energy. 

Acknowledgement: I thank Juliane Behrend for valuable discussions and for her 
interest in this work. 
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